Implementation of quantum algorithms with resonant interactions 
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We propose a scheme for implementing quantum algorithms with resonant interactions. Our scheme only 
requires resonant interactions between two àtoms and a cavity mode, which is simple and feasible. Moreover, 
the implementation would be an important step towards the fabrication of quantum computers in cavity QED 
system. 
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Construction of quantum computer is an enormously ap- 
pealing task because of quantum computational potential to 
perform superfast quantum algorithms. Two classes of quan- 
tum logarithms illustrate the great theoretical promise of 
quantum computers. One is based on Shor's Fourier trans- 
formation including quantum factoring (J], Deutsh-Jozsa log- 
arithm 0] and so on, which are all exponential speedup com- 
pared with the classical algorithms. The other is based on 
Grover's quantum search y|], which is quadratic speedup 
compared with the classical ones. The Grover search algo- 
rithm is very important because many techniques based on 
the search algorithm are universally used in our lives. The 
Grover search algorithm is efficient to look for one item in an 
unsorted database of size TV = 2™ fHIU. Classically, in or- 
der to achieve the task, one needs O(N) queries. However, 
one needs 0(y/~N) queries by the Grover search algorithm. 
Furthermore, the efficiency of the algorithm has been manipu- 
lated experimentally in few-qubit cases via Nuclear Magnètic 
Resonance (NMR) HH] and òptics (JH]. 

Here, we first review the general Grover search algorithm. 
The circuit diagram of the Grover search algorithm with n 
data qubits and one auxiliary working qubit, which can be 
used to search one item from 2" items, is shown in Fig. Q] 
The process can be concluded as the following seven steps: 



1. Prepare the n + 1 qubits, which are in |0)® n |l 



n+l> 



2. Perform the n + 1 Hadamard transformations on the 
n + 1 qubits ; 

3. Apply the oracle. The auxiliary working qubit can be 
omitted after the step; 

4. Perform the n Hadamard transformations on the n data 
qubits; 

5. Apply a phase shift to the data qubits except|0)® n , 
which can be described by the unitary operator 
2|0) lSm (0| — I where / is the identity operation on the 
data qubits; 
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FIG. 1 : The circuit diagram for the Grover search algorithm. H® n 
denotes n Hadamard transformations on the n data qubits. An aux- 
iliary working has been omitted. If x is the one to be searched, 
f(x) = 1, otherwise, f(x) = 0. 



6. Perform the n Hadamard transformations on the n data 
qubits again; 

7. Repeat steps 3^6 with a finite number of times, then 
measure the n data qubits. 

The number of repetitions |@] for obtaining a finite item is 



R = C7(; 



which is bounded above by 
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On the other hand, in the realm of atom, cavity QED tech- 
niques, where àtoms interact with a quantized electromagnètic 
field, have been proved to be a promising candidate for re- 
alizing the quantum processors. Recently, many schemes of 
quantum algorithms have been proposed based on cavity QED 
techniques. For example, Rauschenbeutel et al ifïoll have been 
realized a two-qubit phase gate experimentally with resonant 
interaction of a two-level atom with a cavity mode and Zheng 
ifTTll has realized a two-qubit controlled-phase gate with reso- 
nant interaction of two three-level àtoms with a cavity mode. 
The Deutsh-Jozsa (D-J) logarithm ÍT2I1 and the Grover search 
algorithm lfÏ3Í [l·lll have been realized in cavity QED, and so 
on. 

In this paper, we first mainly propose a simple scheme for 
implementing the Grover search algorithm in cavity QED. 
Comparing the Refs ÍT3I [Till . they are both based on non- 
resonant interactions, our scheme is based on single resonant 
interactions between àtoms and cavity and does not use the 
cavity mode as the data bus. Thus the current scheme is very 
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FIG 2: The circuit diagram for the two-data-qubit Grover search 
algorithm. H denotes Hadamard transformation. \(f>)i and \<f>)2 are 
two data qubits, |e)3 is an auxiliary working qubit. 



simple and the interaction time is very short, which is impor- 
tant in view of decoherence. More importantly, we strictly 
investigate the case of atòmic spontaneous emission and cav- 
ity decay during the interactions. Then we avoid the effect by 
constructing appropriate unitary transformations. Therefore 
our proposal is more approach to real case and can succeed 
with higher fidelity (over 0.99). Here, we only discuss the 
case of two data qubits, where we can search a finite item from 
four items. The circuit diagram of the Grover search with two 
data qubits and one auxiliary working qubit is shown in Fig. 
12 Three-level àtoms are used in this paper and the relevant 
level structure is shown in Fig. [3] The third level \i) is not 
affected during the atom-cavity resonant interaction. Thus we 
consider the case that two àtoms interact with the single mode 
cavity mode, in the interaction picture, the Hamiltonian of the 
atom-cavity interaction can be expressed as (assuming Ti = 1) 
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FIG. 3: The level structure of the àtoms. \g) is the ground state, 
|e) is the excited state. The cavity mode is resonantly coupled to 
the |e) <-> \g) transition. The third level \i) is not affected by the 
interaction. 



So the total state of the àtoms 1 , 2 and 3 becomes 

I0>i23 = ^7^99)12 + \9i)i2 + \eg)u + \ei) 2 )(\g) 3 - |e) 3 )- 

(5) 

Obviously, we have known the four items ( \gg) 12 , \gi) 12 
> \ e g)i2 > |ez)2 ) are stored in the data qubits before applying 
the Oracle. Without loss of generality, we search the state 
\eg)i2 from the four states. However, for the two-data-qubit 
Grover search algorithm, the oracle has effect on the states to 
be searched. The auxiliary working qubit can be discarded at 
this point. 

Secondly, we send àtoms 1 and 2 through the vacuum cav- 
ity, the evolutions are governed by the Hamiltonian of Eq. (Q}, 



H = gi {aïSï + aS+) + .92 (a f ^ + aS+) 



(1) 



where g\ and g 2 are the coupling strength of the àtoms 1, 2 
with the cavity, respectively. ,s + = |e)(<7|, s~ = |g)(e| and \g) 
is the ground state of the àtoms, |e) is the excited state of the 
àtoms, at, a are the creation and annihilation operators of the 
cavity mode. Assume that the cavity mode is initially prepared 
in the vacuum state |0) c . In order to implement the two-data- 
qubit Grover search algorithm, firstly, we prepare àtoms 1, 2 
and 3 in the state 



|e 5 )i 2 |0) c - ÏL{hg lC os(Et) + ^)\eg) 12 \0) c 
Üj £j g\ 

52 [cos(Et) - 1] \ge) u |0) c - isin(Et) \gg) 12 1 l) c }, (6a) 



|ei)i 2 |0) c -» [cos(5ií)|e)i|0) c - isin(git)\g)i\l) c ]\i)2, 

(6b) 



l5fl>i 2 |0) c - \gg) 12 \0)c, 



(6c) 
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(2) 



and send àtoms 1 and 3 through a classical field and choose 
appropriately phase and amplitude, respectively 
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(3a) 



(3b) 



Then we send atom 2 through two classical fields and choose 
appropriately phases and amplitudes 
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where E = \/ g\ + g\. If we choose 

t = —,.92 = V3gi, 
5i 



(6d) 



(7) 



which can be achieved by choosing coupling strengths and 
interaction time appropriately. Thus, we have 
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(8b) 



(8c) 
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|ffï>12|0) c - |fl»>12|0) c . 



(8d) 



Then send atom 2 through two classical fields tuned to the 
transition 



-> |e) 2 , |5>2 < — ► |e) 2 - 
These lead the state of àtoms 1 and 2 to 



(9) 



10) 12 = \{\ag)i2 + \ge)i2 - \eg)i2 + |ee)i 2 ). (10) 

Thirdly, we send àtoms 1 and 2 through a classical field, 
respectively. Choosing appropriately phase and amplitude, let 



\g)i^^(\g)i + \e}i), (i = 1,2) (lla) 



|e) 1 ·-.-^(| ff ) i -|e) i ).(i = l,2) (11b) 

Then we perform a single-qubit operation on atom 2 again 

|e) a -» |i)a- (12) 

Thus Eq. (ÍTÜl l becomes 

10) 12 = 2(b5)i2 + |e.g)i2 - \gi)i2 + \ei)u)- (13) 

In order to achieve the next step (phase transformation), 
we can perform single-qubit operations and controlled-phase 
transformations on the two àtoms as in Eq. (©, which can lead 
Eq. Q3) to 

|0)i2 = ^(\gg)i2 + \ge)i2 - \eg)i 2 - |ee)i 2 ). (14) 

Finally, we perform single-qubit operations on àtoms 1 and 
2 as in Eq. (fTTT i. Thus we obtain the state of àtoms 1 and 2 
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(15) 



We can measure the state of àtoms 1 and 2 by detectors. Ob- 
viously, the state of àtoms 1 and 2 is the result that we want 
to search. If we want to search other states (\gg) 12, \gi)i2 or 
\ei) 12), the main process is the same as above (Eq. ([8]) and Eq. 
(fTTTi). except for some single-qubit operations. 

But in the real processing of resonant interactions, the cav- 
ity decay and the atòmic spontaneous emission are unavoid- 
able. Taking them into consideration, if we choose appropri- 
ate parameters 172 = v3ffi, t = and k = t = 0.1<?i, 
where k is the cavity decay rate and t is the atòmic sponta- 
neous emission rate, the evolution ifTTll of system is similar to 
Eq. ®. While ÍCT^ 20 is added to the \eg) 12 and -|eí)i 2 
compared with the ideal case. The state of Eq. ( ÍT4l i becomes 



;i2 = 



\J\ + 2 x ICTtó + 10"f 
+ 10 _ ^|e5)i2 - 10~TC| ee } 12 ). 



(I35)i2-10 ™\ge)xi 
(16) 



Then we perform the single-qubit operations 
1 



and 
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(10-*/ ao | fl )i + |e)i), (17a) 
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(18b) 



on àtoms 1 and 2, respectively. These lead the state of àtoms 
1 and 2 to \eg)i2, i.e., we can search the state perfectly ( the 
successful possibility and fidelity are both equal to 1.0 ). 

Out of question, the D-J algorithm can be also implemented 
with resonant interactions. The D-J algorithm can distinguish 
the function f(x) between constant and balanced |2]. The 
vàlues of the function f(x) are either or 1 for all possi- 
ble inputs. The vàlues of balance function are equal to 1 for 
half of all the possible inputs, and for the other half. The 
constant is always 1 or for all inputs. Classically, if we 
want to unambiguously distinguish between constant and bal- 
anced function on 2" inputs, we will need 2"/2 + 1 queries 
to achieve the task. While for the D-J algorithm, we will 
need only one query. Here we discuss the two-qubit D-J al- 
gorithm. The state of query and auxiliary working qubit is 
prepared in (|0)j + |1}í)(|0)„ - |l) )/2. Then mapping a 
unitary transformation Uf on the system, which will lead the 
initially state to [(-1)^10), + (-l)« 1 ï|l> i ](|0>„ - |l)„)/2. 
There are four possible transformations to the Uf: (1) for Ufi, 
/(0) = /(l) = 0; (2) for U f2 , /(0) = /(l) = 1; (3) for U f3 , 
/(0) = and /(l) - 1; (4) for U f4 , /(0) = 1 and /(l) = 0. 
After a Hadamard transformation on the query qubit, the state 
of query qubit becomes |/(0) © /(l))- If the function f(x) is 
constant, the state of query qubit becomes |0)j. Otherwise it 
becomes |l)j. 

Obviously implementation of the unitary transformation Uf 
is the key. We prepare two àtoms in 



k)i2 = (|5)i + |e} 1 )(| 5 )2-|e) 2 )/2. 



(19) 



In the case of Uf±, we take no operation on the two àtoms. 
In the case of C//2, we perform a single-qubit rotation on 
atom 2, 



1.9). 



eh. 



(20) 



In the case of Uf3, we first perform a single-qubit rotation 
on atom 2 
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(21) 



Secondly we send àtoms 1 and 2 through a vacuum cavity. 
we can obtain the evolution of Eq. dH), which is governed 
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by the Hamiltonian of Eq. (Q}. Thirdly we perform single- 
qubit rotation on atom 2 of Eq. (|2TT > and another single-qubit 
rotation of Eq. ( f20b . Then we send àtoms 1 and 2 through the 
vacuum cavity again as Eq. dHJ. Finally, we perform single- 
qubit rotation on atom 2 of Eq. (|2TT i. 

In the case of Uf 4, we can achieve the task by the process 
of the case of Uf$ and a single-qubit operation on atom 2. 

Now, we have completed the unitary Uf. Thus the two- 
qubit D-J algorithm will be implemented simply. Moreover, 
the scheme can be generalized to multi-qubit case. 

For the real processing (with cavity decay and the atòmic 
spontaneous emission), we can also achieve the task with suc- 
cessful possibility and fidelity are both equal to 1.0 by choos- 
ing appropriate sing-qubit operations as the implementation of 
Grover search algorithm. This is not only a useful character in 
experimental manipulation but also important for constructing 
real quantum computer. 

Discussion on the feasibility of the current scheme is nec- 
essary. The scheme requires two àtoms in a vacuum cavity 

have different coupling strengths with the cavity mode. The 

2 / 2 

coupling depends on the atòmic positions: g = Vte~ r /w , 
where Í7 is the coupling strength at the cavity center, u is the 
waist of the cavity mode, and r is the distance between the 
atom and the cavity center ÍÏ5ll . The condition g 2 = V3gi 
in our scheme can be satisfied by locating one atom at the 
center of the cavity and locating the other one at the position 
r = a>/n 1 / 2 v / 3. According to the recent experiments with 
Cs àtoms trapped in an optical cavity líïóTl . the condition can 
be obtained. For the resonant cavity, in order to implement 
quantum algorithms successfully, the relationship between 
the interaction time and the excited atom lifetime should be 
taken into consideration. The interaction time should be much 
shorter than that of atom radiation. Hence, atom with a suffi- 
ciently long excited lifetime should be chosen. For Rydberg 
àtoms with principal quantum numbers 50 and 51, the radia- 
tive time is Xi ~ 3 x 10~ 2 s . From the analysis in Ref fÏ7Íl . 
the interaction time is on the order T ~ 2 x 10~ 4 s , which 
is much shorter than the atòmic radiative time. So the con- 
dition can be satisfied by choosing Rydberg àtoms. Further- 
more our scheme requires that two àtoms be simultaneously 
sent through a cavity, otherwise there will be an error. As- 
sume that during the interaction between the two àtoms and 



the cavity, one atom enters the cavity O.Olí sooner than an- 
other atom, with t being the time of each atom staying in the 
cavity. We can obtain the fidelity F ~ 0.999 for Eq. ( fTOt and 
the total fidelity is about 0.998 for the two data qubits grover 
algorithm. Obviously in this case the operation is only slightly 
affected in current schemes. 

Next, one needs to reach the Lamb-Dicke regime in order 
to implement the quantum algorithms. For the state of Eq. 
(O, in the Lamb-Dicke regime, the infidelity caused by the 
spatial extension of the atòmic wave function is about A ~ 
(ka) 2 Tr, where k is the wave vector of the cavity mode and 
a is the spread of the atòmic wave function. Setting A ~ 
0.01, so we have a ~ 0.01A, where A is the wavelength of the 
cavity mode. If the atom trajectòries cross the cavity with the 
deviation of less 0.1 degree from its pre-determined direction, 
we can ensure the fidelity is about 0.999 for Eq. ( fTOb and the 
unitary transformation U / . While in order to maintain g 2 = 

\/3gi in the process of atòmic motion in the cavity, we can 

2 

choose the parameter of cavity z < 0.5z , where z = ^4— 
and 2z is the length of the cavity. We can obtain the error is 
only about 10 ~ 3 . Therefore our scheme is feasible with the 
current cavity QED technology. 

In conclusion, we have proposed the scheme for imple- 
menting the quantum algorithms in cavity QED. Our scheme 
only requires resonant interactions between two àtoms and a 
cavity mode. The interaction time is very short, which is very 
important in view of decoherence. Mean while, even if we 
take the cavity decay and atòmic spontaneous emission into 
consideration, we can still achieve the task perfectly. More- 
over, the implementation of the algorithms would be an impor- 
tant step to scale more complex quantum algorithms and our 
scheme would be very important for constructing real quan- 
tum computer. 
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